Analytic derivation of the map of null rays passing near a naked singularity 
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Recently the energy emission from a naked singularity forming in spherical dust collapse has been 
investigated. This radiation is due to the particle creation in a curved spacetime. In this discussion, 
the central role is played by the mapping formula between the incoming and the outgoing null 
coordinates. For the self-similar model, this mapping formula has been derived analytically. But 
for the model with C°° density profile, the mapping formula has been obtained only numerically. 
In the present paper, we argue that the singular nature of the mapping is determined by the local 
geometry around the point at which the singularity is first formed. If this is the case, it would be 
natural to expect that the mapping formula can be derived analytically. In the present paper, we 
analytically rederive the same mapping formula for the model with C°° density profile that has been 
earlier derived using a numerical technique. 
PACS 04.20.Dw, 04.70.Dy YITP-00-57 



I. INTRODUCTION 

There are many known examples in the literature 
showing that naked singularities can form from regular 
initial data, in the classical gravitational collapse of a 
bounded object. Quantum effects analogous to Hawking 
radiation are expected to play an important role in deter- 
mining the final outcome of such a collapse. A study of 
such effects has received some attention in recent years. 
The most notable feature of these studies is that if a 
massless scalar field is quantized on the background of 
a star forming a naked singularity, the outgoing quan- 
tum flux of the scalar field can be calculated in the ge- 
ometric optics approximation and in 2-d models 
p-p| . This quantum flux can be shown to diverge on the 
Cauchy horizon. An interpretation of this divergence has 
recently been given by || . 

The calculation of the quantum flux in the geomet- 
ric optics approximation crucially relies on the determi- 
nation of the map between ingoing null rays coming in 
from X~ and outgoing null rays going to X + . This map 
can be calculated in a straightforward manner for the 
model of self-similar spherical dust collapse. However, 
when the self-similarity condition is dropped, the calcu- 
lation becomes difficult and the map has so far only been 
determined numerically. 

In this paper we argue that this map is determined, 
not by the global geometry, but by the local geometry 
around the point at which the singularity is first formed. 
It is then possible to determine the map analytically. 

In this paper, we solve the radial null rays in the 
Tolman-Bondi spacetime approximately. These null rays 
naturally defines a map between the incoming and out- 
going null coordinates, u and v. As noted above, this 
map plays a crucial role in calculating the energy emis- 
sion due to the quantum effect in curved spacetime in the 
geometrical optics approximation. 

For the self-similar model with a stellar surface 
on which the background solution is matched to the 
Schwarzschild metric, this map was determined analyti- 
cally 0] . It was shown that the map takes the form of a 
power law : uq — u cx (vq — v ) 7 , where u — uq and v = vq 
are outgoing and incoming null rays passing through the 



point at which the singularity is first formed. 

On the other hand, such a map was not derived analyt- 
ically for the more generic so-called p2-models for which 
we have non-vanishing second derivative of the density 
profile at the center before the formation of naked singu- 
larity. Numerically, the map was shown ||,|| to be given 
as d 2 (u — u)/dv 2 cx (v — v)~ x / 2 . 

In this paper, we would like to show two things. First 
we shall point out that it is not necessary to study the 
global solution for radial null geodesies far from the sin- 
gularity for the purpose of investigating the singular be- 
havior of the map. We shall show that the basic fea- 
ture of the map can be extracted just by considering a 
small region near the point at which the singularity is 
first formed. 

To demonstrate this, we propose an alternative map. 
Let us consider sending radial incoming null rays from 
an observer on a comoving shell. These null rays are 
reflected at the center, and come back to the same co- 
moving observer. A radial null geodesic crosses a co- 
moving shell located at a fixed comoving radius r twice 
before and after the reflection at the center. Thus these 
null rays define a map between the sending time and the 
receiving time measured by the proper time for the co- 
moving observer. We shall show that this map possesses 
the same structure of singularity as the previously inves- 
tigated map between null coordinates. In defining this 
map, we will see that the radius of the comoving shell 
can be chosen arbitrary small. One would expect such 
a result because there is no singular feature in the map 
between ticks on a comoving shell at a finite distance and 
that of the null coordinates naturally defined at infinity. 

Secondly, using the same map, we show that the struc- 
ture of singularity for /?2-models can be obtained analyt- 
ically. 

Throughout this paper, we use the units 8ttG = 1 and 
c= 1. 



II. SELF-SIMILAR CASE 

The marginally bound Tolman-Bondi dust collapse is 
described by the metric 



ds 2 = dt 2 - R' 2 dr 2 - R 2 dQ 2 



(2.1) 



and the evolution of the circumferential radius R(t, r) and 
the density p(r, t) is determined by the Einstein equations 



R' 



F(r) 
R 



F' 
R 2 R' 



(2.2) 



In the above equations, a prime and a dot denote a 
derivative with respect to r and t, respectively. Here 
F(r) is the mass function which in case of the self-similar 
model is given by F(r) — Xr when the scaling is R = r at 
the singular epoch t = 0. Collapse is assumed to begin 
at some epoch t < 0. 

This self-similar model was discussed in Ref. The 
evolution of the circumferential radius is given by 



R 3/2 = r 3/2 f j_ _ <± 

r 



where we have introduced a = §vA. In this model, the 
initial density profile at t = ti„ < near the center is of 
the form 



P = Po + PsR 3 + PeR 6 



(2.4) 



Here the coefficient p% = — (16/3)a~ 3 (— ti n ) 5 is negative, 
in which case the central singularity is known to be glob- 
ally naked for a range of values of A. The values of higher 
order coefficients such as pe and pg are accordingly tuned 
so as to realize a self-similar solution. 

The equation which determines the incoming (upper 
sign) and outgoing (lower sign) radial null geodesies on 
this background geometry is given by 



* = o/ = 1 - \ at l r 

+ dr ' (l-a*/r)V3- 
Following Ref. 0], we introduce the variable 
y = (1 - at/rf 3 . 



Then the above equation is rewritten as 
9y 3 



dr 
r dy 



3y 4 T ay 3 



3y T 2a 



E 



3A1 



(2.5) 



(2.6) 



(2.7) 



where af are the roots of 3y 4 =F ay 3 — 3y =p 2a = 0, and 
the last equality defines the coefficients Af , which satisfy 

Ei=i = 1- Equation ( p.7[ ) is integrated to obtain the 
solution, 



£-n<»-<*> 



±^-3A a 



(2.8) 



Here ro± is the integration constant which parameterizes 
different null rays. For the limit of small r while keeping 
t finite, y goes to oo. Then, in this limit, we find 



r 



1 



1 - af/r 



r 



(2.9) 



Here we have defined £ = — t| r= o- From this relation, 
we can read that tq± = a£. We consider a pair of an 
outgoing radial null ray and an incoming one such that 
the latter is the reflection of the former at the center. 
Then both rays have the same value of £, and hence we 
conclude that 



(2.10) 



For a range of values of A the singularity forming in 
collapse is naked. In this case, the outgoing null ray 
emanating from the point at which the singularity first 
occurs, i.e., r = 0, t = 0, forms the Cauchy horizon. This 
null ray is given by y = a~l where we assume is the 
largest real root among . Similarly, the incoming null 
ray terminating at the first singular point is also given 
by y — where, in the same way, is the largest real 
root among . Since we are considering the null rays 
(2.3) close to the above limiting one, we expand the solution 



around y 



Then we find 



r 



— « (y - of) 



i=l 

*C ± (t ± (r)-t)- 3A i, 
where t±(r) = (r/a)[l — (af) 3 ] and 



3A a 



(2.11) 



^4 



3r(a±) 2 



-3A 4 



IK* 



±\-3A a 



As we have mentioned in the Introduction, we are consid- 
ering a map between the sending time and the receiving 
time measured by an observer on a comoving shell. We 
denote the f orme r time as 1 1 and the latter time as t%. 
Then, from fl2.1CQ and fl2.1l|), we obtain 



{t-(r sh )-t 2 ) 



1/(3A" 



{t+(r sh )-h) A t^, (2.12) 



where r s h is the coordinate radius of the comoving shell 
on which the observer resides. Note that the result does 
not change even if we choose a very small value of r s h- 
By a simple computation, we can show that Af /A± here 
is the quantity denoted by 7 in Ref. Hence, we find 
that the basic property of the map around the singularity 
is maintained in this new calculation, which does not use 
any information about the geometry away from the first 
singular point. 



III. P2-MODEL 

In the preceding section, we considered the self-similar 
model in which the second derivative of the density pro- 
file at the center vanishes, and the third derivative is in a 
certain negative range. In this section, we would like to 
consider more general cases in which the second deriva- 
tive does not vanish. 

What we would like to do here is to deduce an anal- 
ogous map for the model with the_(7°° initial density 
profile considered by Harada et al. 
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p(R) = p + p 2 R 2 + ■■■ 
We write the solution of Eq. (Ej) as 



Ri = l F {r){t-t a {r))\ 



(3.1) 



(3.2) 



where to(r) is an arbitrary function. By using the re- 
maining gauge degrees of freedom we set 



*o(r) 



(3.3) 



Then, the density profile near the center at some epoch 
t = ti n < before hitting the singularity which first 
appears at t = 0, r = 0, is given by 



2 

1 + r 



For small r, we assume F (r) 



F 
F 1 

- ar H 



R'- 



r M for small r at t = i; 



(3.4) 

O(r^). Then 
Req uirin g that the 



second term in the round brackets of Eq.(3.4), r + F/F', 
is proportional to R 2 , we find that p should be chosen 
as 3/2. Then the condition p 2 < is automatic, and the 
singularity at the center turns out to be at least locally 
naked. 

Now let us consider a model truncated at the leading 
order for the expansion of F(r), i.e., we suppose F(r) = 
ar 3/2 Then th e equation which determines the incoming 
and outgoing null geodesies is given by 



dt 
dr 



=F-r 
^6 



1/6. 



1 



ft/r 



(1-t/r) 



1/3 



(3.5) 



where, for convenience, we have performed the following 
rescaling of variables, 



9a 
4 



For small r, the right hand side of the above equa- 
tion goes to zero faster than r . Hence the integration 
constant dominates the solution for small r unless it van- 
ishes. Thus, regarding the right hand side as small, we 
consider an expansion of the solution in the form 



t(r)=£-t 1 (r)-t 2 (r) 



(3.6) 



where £ = — t| r= o as before. 

To keep the notation simple, we introduce the function 



1 - (gf/7) 

(1-S)V3 



(3.7) 



Then, substituting the above ansatz (3.6) into Eq.(3.5) 
we get the equations which determine t(r) order by order 

if = T 7 r ^n-tM, 

dr 6 r 



^6 
7 

^6 



dr 



(3.8) 



Recall that 



W(-{/r) = ^0 
6 czr 



(3.9) 



which relation is not at all accidental but follows from 
the starting expression for the geodesic equation dt/dr = 
^fR'. Using this relation, we can calculate t±(r) as 



h(r) = Tj R(-^,r) = T r 7/6 (l + eA) 2/3 . 
The next order term is evaluated as 



(3.10) 



h = 7 - f dr'r' x l\l + t/r'Y'*F'{-t/r') 
6 Jo 

= 7 -e /6 ^^M 1/3 (l + l/M) 2/3 ^'(-l/M). (3.11) 

Expanding the integrand for large r/£, the asymptotic 
form of the integration is evaluated as 



r/t 



dv^il + l/pf^T'i-l/p) 

C2 _(rM!( 1 + | (OT + 



(3.12) 



Here Ci is the integration constant, which cannot be de- 
termined by integrating the expression expanded for large 
r/£. Except for the term containing C2, all other terms 
are completely determined by the asymptotic expansion 
of the integrand. One may notice that we have not intro- 
duced the corresponding constant for t\{r). In the case 
of t\(r) we can see that this constant vanishes by look- 
ing at the explicit expression which is written in terms 
of R. It will be worth stressing that this cancellation of 
the integration constant is not due to the assumed form 
of R. 

Also for ts(r), we can do a similar calculation as 
T Z^/6 f'* dviA'4 2-T\-\j v ) 
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x / dp^l^l + l/pf^T'i-l/p) 
Jo 

+^J"(-1/,)(1 + 1/,) 4 / 3 
T^ 9 / 6 ±^C 2 ^/V/ 6 (1 + ---) 

Ur 9/6 

T ^62- (1 + -" ) - 



(3.13) 



Besides the terms containing the integration constants 
Cfe, terms in tk have the structure of £ m r6 ~ m+1 (m > 0). 

The region of validity of the above expansion will be 
determined by comparing the first and second term. Re- 
quiring that |ti|/£ <C 1, we find that the condition for the 
validity of the above expansion is satisfied for r <§; £ 6 / 7 . 

For sufficiently small £, within this region of validity, 
there is an overlapping region at which the condition 
£/r <C 1 also holds. When this condition £/r <C 1 is 
satisfied, we can expand the geodesic equation as 



dt_ = 7 1/6 / _ 2t_ 

dr 6 V 21r 63r 2 



(3.14) 
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Then, it can be integrated iteratively as 
t-D= = r 7 / 6 




where is the integration constant. Terms in the 
(I + l)-th line in the right hand side have the struc- 
ture of D l -^ n r 1 e - ~ l (I > n > 0). Comparing them 
with the terms in tk, we find the correspondence with 
I = m + k — 1, n = k. Corresponding terms in two differ- 
ent expressions have the same coefficient. In other words, 
the first term in each line of the right hand side basically 
corresponds to h(r), the second term to t 2 {r) and so on. 
The relation between and £ is found by comparing 
the coefficient of the r-independent term. Then we have 

d t + o(d 2 t ) = -t + 7 -c 2 e /& t lc 3 e /e 

+ ^ 4 £ 10/6 + O(£ 11/6 ). (3.16) 

Eliminating £ from this equation, we obtain the relation 
between D + and _D_, 

D-=D + + 7 -C 3 (-D + f/ 2 + 0{D^ ,& ). (3.17) 

This mapping formula has the same structure of singular- 
ity as that obtained for the map between incoming and 
outgoing null coordinates by Harada, Iguchi and Nakao 
numerically. 



IV. INTERPRETATION 



In the £ — > limit, the region of validity of the solution 
for small r shrinks to zero. Hence, the solution for small 
r cannot describe the limiting behavior for any value of 
r s h- 

Anyway, the basic feature of the map is derived for ar- 
bitrarily small r s h . This fact proves that it is determined 
by the geometry just around the central singularity as ex- 
pected. This is unlike the case of the black hole for which 
the map is not determined by the local geometry near the 
center. Instead it is determined by the relationship be- 
tween the regular interior coordinates and the external 
coordinates, near the event horizon. 

Here we have considered a specific model for the initial 
density profile. However, we can expect that the basic 
structure of the map will be independent of the details of 
the model because it is determined by the geometry just 
near the point at which the singularity is first formed. 
In this small region, we would be able to neglect the 
higher order terms in the expansion of the density pro- 
file, although further consideration is necessary to give 
a rigorous proof for the statement that the structure of 
singularity in the mapping formula is totally determined 
by the leading term in the expansion of the initial density 
profile for general spherical dust collapse. 

In this paper we have discussed models of spherically 
symmetric dust collapse. However, it would be of interest 
to generalize the prescription described in this paper to 
other matter models and to non-spherical collapse. 
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The map that we have obtained in this paper is 
the map between the ingoing and outgoing intersection 
points of null rays measured by the proper time of the co- 
moving observer. This mapping is not exactly the same 
as the mapping of the incoming and outgoing null coordi- 
nates naturally defined for the observers at null infinities. 
However, as we mentioned in the Introduction, the sin- 
gular behavior of these two mappings is expected to be 
common because there is no singular behavior in the null 
rays connecting the shell at a fixed r and the shell at 
infinity. 

Then, in order to investigate the structure of the sin- 
gular mapping, is it possible to set the location of shell 
r s h arbitrarily small? The answer is basically yes. In the 
above derivation, the only constraint on r s h is r s h ^ 
but this constraint does not prevent us from choosing 
r s h arbitrary small. Since we are interested in the lim- 
iting behavior when the time of reflection t = — £ comes 
very close to the singularity t = 0, we can assume that 
£ is arbitrary small. In the preceding section, we consid- 
ered a matching of two different expressions for the null 
geodesies. One is valid for small r and the other for large 
r. At this moment, one may wonder why it was neces- 
sary to consider this matching although we are allowed 
to choose r s h arbitrarily small. The reason is as follows. 
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